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Recent results on the system size dependence of net-baryon and hyperon produc- 
tion as measured at the CERN SPS are discussed. The observed -/Vpart dependences 

O ■ 

CN I of yields, but also of dynamical properties, such as average transverse momenta, can 

Q ■ be described in the context of the core corona approach. Other observables, such 

n. 

as antiproton yields and net-protons at forward rapidities, do not follow the pre- 

I dictions of this model. Possible implications for a search for a critical point in the 

^ ' QCD phase diagram are discussed. Event-by-event fluctuations of the relative core to 

D ■ 

corona source contributions might influence fluctuation observables (e.g. multiplicity 

O _ 

^ I fluctuations). The magnitude of this effect is investigated. 

(N 
> 

. 

^ . The system size dependence of particle production in heavy ion collisions can be measured 

in two ways: either by comparing central reactions of nuclei of different size, or by analyzing 
centrality selected minimum bias collisions of nuclei of fixed size. It turns out that many 
^ ! observables follow a very similar A^part dependence, if their centrality dependence is studied. 

Usually a rapid change for very peripheral reactions is followed by only a weak variation 
towards very central collisions. This indicates that many features of system size dependences 
are dominated by the same underlying geometrically effects. It is therefore important to find 
a reliable way of modeling them in order to extract non-trivial effects. 



I. INTRODUCTION 
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II. CORE CORONA SEPARATION 



A quite successful way of describing the main features of system size dependences of 



a. 



In this model a nucleus-nucleus 



particle production is the core corona approach 
collision is decomposed into a central core, which corresponds to the large fireball produced 
in central A+A collisions, and a surrounding corona, which is equivalent to independent 
nucleon-nucleon reactions. To quantify the relative contribution of the two components the 
fraction of nucleons that scatter more than once, f(A'part); can be used. /(iVpaj-t) can 
simply be calculated within a Glauber model jsl [g]. This quantity allows for a natural 
interpolation between the yields Y measured in elementary p+p (= Icorona) and in central 
nucleus-nucleus collisions (= Icore): 



y{Npe,rt) = iVpart [/(^Vpart) W + (1 - /(iVpart)) >^corona] (1) 

The left panel of Fig. [1] shows /(A^part) for various symmetric reaction systems. Due to 
their different surface to volume ratio, the N-p^^j-i dependence of /(A^part) is much steeper for 
smaller reaction systems than for larger ones. This feature of the model allows, for instance, 
to describe the centrality dependences of strange particle production measured for Cu+Cu 
and Au+Au collisions at RHIC [7]. An important feature of the curves shown in the left 
panel of Fig. [1] is that the maximum value of /(A^part) depends on the size of the nuclei. 
While for very central Pb+Pb collisions /max ~ 0.9 can be reached, the maximum value for 
very central C+C reactions is significantly lower: /max ~ 0.65. This A dependence of /max 
is illustrated by the dashed red line. It has the important consequence that the different 
relative core corona contributions have to be taken into account when comparing even very 
central nucleus-nucleus collisions of different size, which usually are considered to represent 
only the core-like fireball. This effect might explain any observed system size dependence of 
the chemical freeze-out parameters T and /Ib, as derived from statistical model fits to central 
A+A collisions for nuclei of different size (e.g. 8|). As a consequence of this it therefore 
follows that the system size does not provide a good control parameter to probe different 
regions of the QCD phase diagram via the production of fireballs of different temperature 



Is 



9[ . It is thus more likely that one only observes a change in the relative admixture of central 
fireball (core) and peripheral p+p like corona, whose freeze-out parameters are different, but 
independent of the size of the involved nuclei. In any case the underlying effect from the 
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core corona separation should be taken into account when chemical freeze-out parameters 
are extracted. 

Another interesting aspect of the core corona approach is shown in the right panel of 
Fig. [H Here, the A^part dependence of /(A'part) is shown for asymmetric collisions. In this 
case quite distinct centrality dependences can be observed. While for symmetric collisions 
this dependence is following a continuous rise (left panel of Fig. [1]), for asymmetric collision 
systems with a small projectile nucleus (e.g. 0+Pb or Si+Pb) a rapid rise followed by 
a maximum and a subsequent decrease of /(Apart) is seen. Therefore, a similar centrality 
dependence of particle yields can be expected for these type of colhsions and its measurement 
would constitute a test for the validity of the core-corona model [9|. 



III. STRANGENESS PRODUCTION 



The left panel of Fig. [2] shows the system size dependence of the strangeness enhancement 



factors E measured by the NA49 collaboration at ^/s^ 



17.3 GeV 
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ll|. For this 



measurement E is defined relative to p+p reactions as baseline: 



E 



dAr(Pb+Pb) 



V(iVpart) 



dy 



y=0 




1 dAr(p+p) 
2 



dy 



(2) 



As a general feature a rapid rise of E for A^part < 60 is seen, which then turns into a 
slow increase with A^part ^) ^ Q, or even into a saturation (A). Comparable 
observations have been made by the NA57 experiment I4] at the same y/s^ and by the 
STAR collaboration at y/s^ = 200 GeV 13|]. 

In the context of statistical models, the system size dependence of strange particle pro- 
duction was expected to be described by the transition from a canonical to a grand canonical 
ensemble [l^ . This effect provides a natural explanation of the hierarchy of the suppression 
pattern, which depends on the strangeness content of the particles {E{A) < E(E) < E{Q)). 
However, the onset of the enhancement is predicted to happen for smaller systems (A^part ~ 
30), than observed in the data. Also, this model does not explain the slow rise with Apart 
seen for S and Q for mid-central and central collisions. However, within the core corona 
model this behavior can be described naturally, as described in the previous section. This 
is illustrated by the solid lines shown in the left panel of Fig. |2l 

It is interesting to note that this approach not only matches the observed centrality 
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dependences of yields, but does also describe dynamical quantities such as the average 
transverse momenta (m^) — mn - This is shown in the right panel of Fig. |2] for strange 
particles and (anti-)protons jlo[ fisj. A similar observation has been made for elliptic flow 



at RHIC 



IV. CENTRALITY DEPENDENCE OF NET PROTONS 



Recent data on proton and antiproton production allow to study the system size depen- 



dence of stopping [17|. In this analysis the rapidity spectra of net-protons (p — p) have 
been determined in several bins of centrality for minimum bias Pb+Pb collisions at 40A 
and 158^4 GeV. Based on this data the centrality evolution of net-proton yields at different 
rapidities can be investigated. The left panel of Fig. |3] shows a comparison of the net-proton 
yields normalized by A^part (1/^part dn/dy{p — p)) as measured at midrapidity and at for- 
ward rapidity at 158 A GeV. While the midrapidity net-proton yields are slightly rising from 
peripheral to central Pb+Pb collisions, the forward yields exhibit a significantly different 
centrality dependence. Here the normalized yields are constant from very peripheral on or 
even drop towards very central collisions. In contrast to the net-protons at midrapidity, the 
forward yields can therefore also not be described by the core corona approach. This implies 
that for the system size dependence of stopping, i.e. the processes that transport baryon 
number along the longitudinal axis and which dominate the net-proton spectra at forward 
rapidities, the distinction between core and corona is not important. These processes do not 
seem to depend to strongly on the number of collisions that a single nucleon experiences. 

For the yield of antiprotons a different centrality dependence is observed than for other 
particles |17] (right panel of Fig. |3]). In this case the yield per participant is decreasing from 
p+p collisions towards central Pb+Pb reactions. This is most likely caused by absorption 
antibaryons in the fireball medium. Also this effect is difficult to describe in the core corona 
framework, since it rather depends on the average path length of the antiprotons inside the 
fireball medium than on relative size of core and corona. The resulting antiproton yields will 
thus be rather a complicated interplay between the system size dependences of production 
and absorption mechanism. 
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V. MULTIPLICITY FLUCTUATIONS 



The system size dependence of multiplicity fluctuations was measured by the NA49 ex- 



periment [l9l]. The fluctuations, quantifled by the variance of the event-by-event measured 
multiplicity n normalized by its mean u = Va.T{n)/{n), exhibit a quite distinct centrality 
dependence. While at y/s^ = 17.3 GeV w ~ 1 is found for p+p collisions, very peripheral 
Pb+Pb reactions result in much higher fluctuations (w ~ 3 for all charged particles) which 
then decrease towards very central collisions where u close to unity is observed again. In 
the following we investigate to what extend these fluctuations could be caused by geometry 
fluctuations that are intrinsic to the core corona approach. This source of fluctuations is of 
special interest, since it will be present even in the idealized situation of a perfect centrality 
selection, corresponding to the hypothetical situation where the number of participants could 
be flxed exactly in the experiment. Also in the case iVp^rt = flxed, the relative contribution 
of the core and the corona, as given by /(iVpart); "wiH still fluctuate from event to event. 
The upper left panel of Fig. H] shows the distribution of / versus iVpart^ calculated within 
a Glauber model js, The event-by-event distribution of / is much wider for peripheral 
reactions than for central ones. The resulting variance Var(/) (see upper right panel of 
Fig. H]) has a A^part dependence that is quite similar to the one observed for multiplicity 
fluctuations. Since the particle yield per participating nucleon is smaller for the corona than 
for the core part, the fluctuations in / will translate into multiplicity fluctuations that should 
be part of the observed ones. To evaluate how large this contribution can be, a simple Monte 
Carlo model was constructed. In this model the multiplicity of a given event is calculated 
according to Eq. ([T]). The event-by-event yields of the core and corona sources, Fcore and 
^corona? are generated from Poisson distributions, whose means have been adjusted to the 
measured midrapidity yields per participating nucleon of charged pions in central Pb-|-Pb 



(core) 20| and p+p (corona) 2l| collisions. The value of / for a given A'^part is sampled from 
the distribution shown in Fig. HI In the lower left panel of Fig. Hlthe result of this simulation 
is shown. This way one obtains multiplicity fluctuations whose centrality dependence has 
the same shape as the one observed in the data. However, the magnitude of the fluctuations 
is much higher in the data than what results from this model study (see the comparison in 
the lower right panel of Fig. H]). The difference between the yields Icore and Icorona would 
need to be two orders of magnitude higher to reach the level of the measurement. 
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VI. CONCLUSIONS 



Many observables measured in heavy ion reactions turn out to follow a very similar 
centrality dependence, which is characterized by a rapid rise in the region of very peripheral 
coUisions followed by a significantly slower increase towards very central ones. This include 
the yields per participant of strange particles, but also their mean transverse mass and even 
elliptic flow V2- It is found that the core corona model provides a very effective way to 
describe the main features of this geometry driven behavior and thus constitutes a good 
comparison baseline which allows to separate non-trivial effects. Observed exceptions to 
this general picture are, for instance, the centrality dependence of antiproton yields and net- 
proton dn/dy at forward rapidities, where more complicated physical processes are involved 
(absorption of antibaryons, baryon number transport) which cannot easily be described 
by the core corona approach. For asymmetric collisions a signiflcantly different centrality 
dependence than for symmetric ones is predicted by this model, which is characterized by a 
maximum for mid-central and a drop towards central collisions. Therefore the investigation 
of asymmetric collisions might provide a good test for the core corona approach. Another 
consequence of this model is that the observed change of the chemical freeze-out temperature 
as derived from statistical model analyses is basically only due to a change of the relative 
admixture of the core and corona contributions. Therefore, the system size might not provide 
a good control parameter to scan the QCD phase diagram, since the properties of the central 
core fireball might in fact not change and the system always follows the same trajectories 
through the phase diagram. Event-by-event fluctuations of the relative core and corona 
contribution, which must be present even in the case of an ideal centrality selection (i.e. 
A^part = flxed), do contribute to multiplicity fluctuations. However, the magnitude of this 
effect is way too small to explain the measured data. Still, further studies might help to 
understand the system size dependences of fluctuations of non-intensive quantities, such as 
K/tt fluctuations. 
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FIGURE CAPTIONS 

FIG. 1: The fraction /(A^part) of participating nucleons that scatter more than once as a func- 
tion of the number of participants A^part- /(^part) calculated within a Glauber 
model 5|, l6|. The left panel shows results for symmetric systems, the right panel for 
asymmetric ones. The dashed lines connect the values for the most central collisions. 

FIG. 2: Left: The midrapidity yields of A, A, S~, and f2~ + r2+ per wounded nucleon relative to 
)+p yields for central C+C, Si+Si and minimum bias Pb+Pb reactions at 158^4 GeV 
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Uj. Right: The {rrif^) — niQ values of strange particles and (anti-)protons at mid- 
rapidity for Pb+Pb collisions at 40^4 and 158^4 GeV, as well as for near-central C+C 
and Si+Si reactions at 158^4 GeV 10|. The (anti-)proton data are taken from 15 |. 
Also shown are the results from a fit for A and A with the core-corona approach (solid 
lines) . 

FIG. 3: Left: The net-proton yields per number of wounded nucleons at forward rapidity and 
at midrapidity for centrality selected Pb+Pb collisions at 158A GeV. Right: The total 
antiproton and antilambda yields per number of wounded nucleon at midrapidity for 



centrahty selected Pb+Pb collisions at 158 A GeV [l7| 



FIG. 4: The distribution of the fraction of nucleons that scatter more than once /(Apart) 
(upper left) and the corresponding variance Var(/) (upper right), both as a function 
of A^part calculated for Pb+Pb collisions at 158A GeV. The lower left panel shows 
the multiplicity fluctuations generated with the MC model (see text for details). The 
red circles display the result if event-by-event fluctuations of / are included, while 
the blue boxes correspond to the case where they are switched off. The lower right 
panel shows multiplicity fluctuations for all charges versus^ the number of projectile 
participants as measured in Pb+Pb collisions at 158 A GeV [19] (blue boxes). The red 
line is the same MC result as shown in the left panel. 



